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Spin Glass. 

C. Amoruso and A.K. Hartmann 

Institut fur Theoretische Physik, Universitat Gottingen, 
Tammannstrasse 1, 37077 Gottingen, Germany 

M.A. Moore 

School of Physics and Astronomy, University of Manchester, Manchester M13 9PL, U.K. 

Energy barriers determine the dynamics in many physical systems like structural glasses, disor- 
dered spin systems or proteins. Here we present an approach, which is based on subdividing the 
configuration space in a hierarchical manner, leading to upper and lower bounds for the energy bar- 
rier separating two configurations. The fundamental operation is to perform a constrained energy 
optimization, where the degree of constraintness increases with the level in the hierarchy. 

As application, we consider Ising spin glasses, where the energy barrier which needs to be sur- 
mounted in order to flip a compact region of spins of linear dimension L are expected to scale as 
V* '. The exponent ip is very hard to estimate from experimental and simulation studies. By using 
the new approach, applying efficient combinatorial matching algorithms, we are able to give the first 
non-trivial numerical bounds 0.25 < ip < 0.54 for the two-dimensional Ising spin glass. 

PACS numbers: 75.10.Nr, 75.40.Mg, 02.10.Jf 
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A unifying concept in the physics of disordered systems 
is the notion of energy landscapes. The dynamics of these 
systems is determined by (free) energy barriers. Promi- 
nent examples include spin glasses [lj, structural glasses 



and folding of proteins [3j. The barrier problem is 
also related to theoretical computer science, because it 
belongs 4] to the fundamental class of nondeterministic 
polynomial (NP) hard problems ||, hence all known al- 
gorithms determining lowest barriers take a time growing 
exponentially with system size. 

In this paper, we will present a hierarchical approach 
to calculate minimum barriers, which is based on the 
application of combinatorial optimization algorithms [6| . 
As application, we apply it to the prototypical two- 
dimensional Ising spin glass, where the value of the 
barrier-height exponent could not be determined so far. 

In the droplet Q or scaling theory of Ising spin glasses 
the low-energy excitations are compact droplets. The 
creation of a droplet (a region of reversed spins) results in 
the formation of a domain wall around the reversed spins 
and the energy of the droplet scales as L e where L is the 
linear extension of the droplet. The system orders at low 
temperatures only if 9 > 0. The domain wall is fractal 
and has size L ds , with d — 1 < ds < d. The dynamics 
of the system are controlled by the height of the barriers 
which have to be crossed to create such droplets. It is 
generally assumed that the barrier to be surmounted to 
create a droplet of linear extent L has an energy which 
scales as where 6 < ip < d — 1 for dimension d. The 
argument for these inequalities goes as follow: the barrier 
must be at least as large as the energy required to create 
the domain, hence 9 < ip; the upper limit is due to the 
fact that the barrier must be lower than the energy of 



the same number of reversed spins. While for a directed 
polymer in a random system it was shown that ip — 9, 
for Ising spin glasses ip appears to be an independent 
exponent. 

While many numerical studies have been done to cal- 
culate the exponent 9, we are aware of only one direct 
numerical estimate of ip in two dimension It was 
equal to the upper bound i.e, ip = d — 1, but only small 
systems were studied {L < 6). As mentioned, the bar- 
rier problem is NP hard and it has been difficult even to 
find good approximate algorithms for barriers. This is 
presumably the reason for the paucity of studies of this 
exponent. Recently Drossel and Moore instead of at- 
tempting to calculate the barrier exactly, placed bounds 
on its energy They showed that for the hierarchi- 
cal (i.e. Berker) lattice there exists a lower bound on ip 
which is the same as the upper bound d — 1 and hence 
they concluded that ip = d— 1. In the same spirit as their 
work, we will also obtain upper and lower bounds on the 
energy of the barriers but for the more physically relevant 
case of the square lattice. Rather than study droplets, 
we will consider the computationally simpler (but equiva- 
lent) task of determining upper and lower bounds for the 
barrier separating the two ground states (GSs) related 
by a flip of all spins. The algorithm is based on subdi- 
viding the configuration space in a hierarchical manner 
and performing constrained energy optimizations, where 
the degree of constraintness increases with the level in 
the hierarchy. This approach can be applied in general 
to obtain minimum barrier for many problems, but for 
the convenience of the reader, we formulate it here for 
the case of the two-dimensional Ising spin glass. 

The model consists of N — L 2 spins Si = ±1 on a 
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direction and free boundary conditions in the y direction. 
The Hamiltonian is 

H = -^J ij SiS j , (1) 

m 

where the sum runs over all pairs of nearest neighbors (ij) 
and the are the quenched random spin couplings. We 
will consider a Gaussian distribution of couplings with 
zero mean and unit width. We define the barrier to be 
the energy required to invert all spins in the GS within 
the context of single spin flip dynamics, as in y|. Each 
trajectory of reversal is characterized by the highest max- 
imum in the trajectory having energy SE max , in excess 
of the GS energy, so that Eb = min(5 E max ) , where the 
minimum has to be computed over all the possible tra- 
jectories. 

Thus the first task is to compute the GS {a } of 
the spin glass; in general already this problem is NP 
hard, but for planar graphs it can be mapped on to the 
minimum-weight perfect matching problem, which is solv- 
able in polynomial time. The algorithms for minimum- 
weight perfect matchings [lfl Hll | are among the most 
complicated algorithms for polynomial problems but the 
LEDA library offers a very efficient implementation ^3] • 
For the details of the method, see Refs. @, 

We will now show how is it possible to extract bounds 
on Eb- In general, a lower bound on the energy barriers 
can be obtained, if one considers a sequence of constraints 
fixing some degrees of freedeom leading from one config- 
uration to the other and calculating the GS under each 
constraint. Here, we force a domain wall, which seperates 
spins having orientation {a } from spin having orienta- 
tion {—a }, through constraining successive bonds along 
the x-axis. Thus, for each realization of the bonds Jij, 
a GS is first calculated with free boundary conditions in 
both directions. Then "hard" bonds |l3( are introduced, 
i.e. bonds with a high value of the absolute strength, 
compatible with the relative GS orientation of the adja- 
cent spins, in a line which runs from the left border to 
the right border. Next the sign of exactly one hard bond 
on this line is inverted and the GS of the new realization 
is calculated. L times exactly one bond is inverted. The 
minimum of these L energies is the domain wall associ- 
ated with the exponent 9. The maximum can be asso- 
ciate with a lower bound on ijj, because one is calculating 
the minimum with the highest energy in the sequence 
from putting the hard bonds at successive places on the 
x-axis. The true barrier height is associated with flip- 
ping one spin at a time and one might have to climb an 
even higher barrier to reach the maximum domain wall 
energy starting from the domain wall of minimum en- 
ergy. We will next explain how to place an upper bound 
on the barrier height before discussing the results of the 
numerical work. 

The basic idea of the general algorithm is to continue 



i.e. one introduces additional constraints between all ad- 
jacent pairs of constrained GSs obtained above. This is 
continued recursively, until the level of changes of single 
degrees is reached. Hence, one obtains a path in config- 
uration space where the maximum energy configuration 
yields an upper bound for the minimum energy barrier. 
I.e. For the spin glass we do this by introducing hard 
bonds along the line corresponding to the maximum spac- 
ing between two "adjacent" domain walls (either in the 
x or the y direction) and by forcing the first domain wall 
through each of these bonds, see Fig. ^ This procedure 
is repeated itcratively until two consecutive domain walls 
differ only by a unique sequence of spin flips. Of course, 
since we are minimizing over a restricted set of paths, 
there could be some lower path over the energy land- 
scape, hence our result will be an upper bound on the 
energy barrier. However, since the upper bound makes 
an attempt to estimate the energy needed to move the 
domain wall to successive places along the x-axis by a 
sequence of single spin flips, we believe the true value of 
V> is likely to be nearer the upper bound than the lower 
bound. 




FIG. 1: Example showing how the domain wall is shifted. In 
(a) the black and the gray area represent the domain walls 
created by changing the sign of two consecutive hard bonds 
in the top (thick) line. The maximum distance between the 
two domain walls is determined and hard bonds are intro- 
duced along this line (thick line at the middle). Then the 
sign of each of these bonds is inverted exactly once (wiggled 
segment) and the new GS is calculated, so that the domain 
wall is forced to pass through them . In (b) the grey area 
now represents the configuration obtained when changing the 
sign of the leftmost hard bond and calculating the new GS. 
Again hard bonds are introduced along the line of maximum 
distance and the new GS (c) is computed for the realization in 
which the leftmost hard bond is inverted, as in (a). Now the 
distance between the two domain walls is at most one spin, 
so this spin is now flipped (yielding one step of the whole 
seauence") and the resulting configuration is shown in (d). 
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For each flipped spin the energy of the configuration 
is calculated, and at the end the maximum among the 
iV energies is considered. In this way we give an explicit 
rule to build a sequence of spin-flip taking one from the 
GS configuration {a } to {—a }. An example of the cor- 
responding energy landscape sampled by this sequence is 
shown in Fig. [21 
One has to be careful to exclude the domain walls which 
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FIG. 2: An example of the energy landscape explored by our 
algorithm for L — 14; the exponent ip is associated with the 
scaling of the maximum energy found in the sequence. 

are not crossing the system and for which the width L 
is not an appropriate measure of their scale. Those con- 
figurations can appear in the first steps of the algorithm 
described above because of the use of free boundary con- 
ditions. Fortunately most of the domain walls do span 
the system. We have restricted our analysis to those con- 
figurations. 
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FIG. 3: Upper line: maximum over L energies in the sequence 
from putting hard bonds at successive places on the x axis 
plotted as a function of the system size. Lower line: maximum 
energy in the sequence of single spin-flip built as explained 
above. 

We applied our method for system sizes in the range 
8 < L < 40, and for each L used 1000 independent re- 
alizations of the bonds Jy. The two quantities AE rn i„ 
(the lower bound) and AE max (the upper bound) as a 




FIG. 4: Rescaled probability distribution of the lower bound 
for the barriers heights with ip = 0.25 



tities can be fitted by an algebraic function ~ Lr , where 



0.25 ± 0.01 and i> up = 0.54 ± 0.01. There could of 



course be a procedure to obtain a better path in configu- 
ration space, yielding a barrier which grows more slowly 
than L 5i . However this is to our knowledge the only nu- 
merical estimate of the energy barriers exponent giving 
an upper bound different from the trivial one d — 1. 

There are many experimental estimate of iJj in the lit- 
erature, and there is disagreement as to its value. Dekker 
et al |l-"il| reported experimental verification of activated 
dynamics in a d — 2 system obtaining ip = 0.9. Schins et 
al. ^(| find ip — 1-0 ± 0.1 by studying aging via the low- 
frequency ac susceptibility giving credence to the claim 
that ip is equal to its upper bound. However, Dupius et 
al [lg extracted 0.3 < -0 < 0.7 in d = 2 and i/i ~ 1.1 
in d = 3. The numerical study of Gawron et al. gave 
ip = 0.9 ± 0.1 in 2 dimension, while that of Berthier and 
Bouchaud |l8| had ip ~ 1.0 and ip ~ 2.3 in 3 and 4 
dimensions respectively. We remark that our optimiza- 
tion algorithms do not suffer from equilibration problems, 
they are exact and they allow the study of large systems 
at least in two dimensions. Besides, by iteratively apply- 
ing our matching algorithm, we are producing large scale 
changes to the domain wall, which is how ip can become 
less than d — 1. (We have studied the barrier associated 
with a straight domain wall and found that the expected 
linear behavior is approached by L ~ 40). 

In Figs.0]and|Hlthe rescaled probability distributions of 
the barriers heights corresponding to both the lower and 
upper bounds which we have found, averaged over the 
disorder, is shown for different sizes L: in both cases it 
seems to approach a fixed shape in the thermodynamical 
limit (although the approach to a scaling collapse is faster 
for the upper bound). 

Notice that both these distributions have long tails in- 
dicating the existence of some very large barriers, which 
implies that some spin configurations could be very long- 
lived. 

We also investigated geometrical properties of the do- 
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FIG. 5: Rescaled probability distribution of the upper bound 
for the barriers heights corresponding to if) = 0.54 

boundaries of the low-energy excitations are fractal with 
a an average length along the perimeter of order P ~ L ds , 
with d — 1 < ds < d. We have computed the fractal 
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FIG. 6: The perimeter (length) of the surface of the domain 
wall corresponding to the upper bound (maximum) energy. 
Data are fitted very well by a power law P ~ L ds with ds = 
1.27 ±0.01. 

dimension of the domain wall ds corresponding to the 
configuration of maximum energy. We have found the 
interesting result that ds = 1-27 ± 0.01, which is the 
same value for domain walls of minimum energy. 

To conclude we have introduced a hierarchical algo- 
rithm to compute upper and lower bounds on energy 
barriers in disordered system. Applying the algorithm 
to two-dimensional Ising spin glasses, we find that the 
minimum barrier energy is bounded above by the scale 
L 54 and below by L 024 . Our numerical upper bound 
is significantly less than the rigorous upper bound (and 
the value for the hierarchical lattice), d — 1. Hence, we 
suspect that ip is an entirely non-trivial exponent. At the 
present time there seems no theoretical approach which 
might explain a non-trivial value of if) 
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